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Abstract 



C ■ We extend an effective Lagrangian embodying broken scale and chiral symme- 



try to include explicit chiral symmetry breaking and an additional chiral invariant 
term which allows for an axial coupling constant greater than unity. We also 
include a chiral Lagrangian for the isotriplet vector mesons which leads to a renor- 
malization of the pion field. The properties of nuclear matter and nuclei, low 
energy vrA^ scattering and the behavior of quantities such as the pion mass and 
axial coupling at finite density are discussed. 



1 Introduction 



Since spontaneously broken chiral symmetry is thought to be a fundamental 
feature of low energy effective Lagrangians, it is remarkable that it is difficult 
to describe nuclear matter and nuclei with such a Lagrangian. For example, 
it has long been known |Q that the standard Lagrangian of the linear sigma 
model, supplemented by the repulsion of the uo meson, does not yield a normal 
saturating equation of state for nuclear matter. In order to describe nuclear 
matter, the Lagrangian can be modified by generating the vector meson 
masses through coupling to the scalar a field or by including one-loop vacuum 
corrections. However, it is found |0, ^, ^ that such models fail to reproduce 
the observed properties of nuclei. 

In a previous paper 0, hereinafter referred to as I, we were able to sat- 
isfactorily describe nuclear matter and finite nuclei with an effective La- 
grangian which incorporated broken scale symmetry in addition to sponta- 
neously broken chiral symmetry, as suggested by quantum chromo dynamics 
(QCD). In order to obtain good phenomenology it was necessary to gen- 
erate the vector meson masses by coupling to the glueball field, 0, and to 
discard the conventional "Mexican hat" potential of the linear sigma model, 

(a^ + TT^ — f^f. In fact by introducing an {uj^uo^^Y term in the Lagrangian 
and adjusting the coupling, we obtained results of a quality comparable to 
those of the standard Walecka model with non-linear and terms [^]. 
In our model the scalar masses are 1.5 GeV and approximately 0.5 GeV for 
the states which are predominantly (j) and o", respectively. Here the a is 
the chiral partner of the pion and also provides intermediate range attrac- 
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tion. This view of the a meson finds support in the work of the Brooklyn 
group using the Nambu-Jona-Lasinio model to study correlated two-pion 
exchange. They find that for the spacelike momentum transfers of interest 
here, the T matrix can be represented by an effective low-mass sigma me- 
son which is the chiral partner of the pion, while in the timelike region a 
physical, low-mass scalar particle is not obtained. On the other hand Furn- 
stahl et al. ^ eliminate the scalar chiral partner of the pion by using the 
non-linear model and consider that a separate low-mass scalar generates the 
intermediate range attraction. They also incorporate broken scale invariance 
and obtain an effective Lagrangian which gives a good account of nuclear 
matter and finite nuclei. (The experimental situation in the scalar sector 
unfortunately remains confused [§].) 

The purpose of the present paper is to further test our Lagrangian by 
explicitly studying pions, since they made no contribution in I. This requires 
some extension of the Lagrangian. Firstly, explicit chiral symmetry breaking 
is needed to endow the pion with a mass. Secondly, in I the axial coupling 
constant, qa-, was unity, rather than the actual value of 1.26; this defect 
can be corrected by introducing an additional chiral invariant term in the 
Lagrangian, as pointed out by Lee 0. Thirdly, if the isotriplet vector and 
axial vector mesons (p and ai) are included in a linear chiral model, along 
with the scalar and pseudoscalar mesons, vr — ai mixing ensues |jlO| and 



this we explicitly take into account. The theoretical development is given 
in Sec. 2, while in Sec. 3 we give results for nuclei, vrA^ scattering and the 
predicted density dependence of quantities such as qa and the pion mass. 
Our conclusions are presented in Sec. 4. 



2 Theory 



We begin by recalling the effective Lagrangian employed in I. It was written 
in the form 

£1 = £0 - , (1) 

where we have schematically separated out the scale and chiral invariant 
part, from the potential Vq which induces the breaking of scale and 
chiral invariance. We write Cq in terms of the nucleon field A^, the chiral- 
invariant combination of sigma and pi fields, a and tt, the glueball field and 
the field of the omega vector meson cj^ which is a chiral singlet. Specifically 

+ [(6-4)^^1' + N[Y{^^, - g^u^) - g{a + ztt ■ t^,)]N . (2) 

Here the field strength tensor is defined in the usual way, cj^,^ = d^u,^ — d^ujfj, 
and we have included a term quartic in the vector meson field (note that since 
the coefficient of this term is positive, the violation of causality discussed in 



Ref. [0 does not arise here). In principle the scale invariant 00 mass term 
could be generated by coupling to the or the cr field, however the latter 

choice does not yield a good phenomenology for nuclei 0, |^, Q so we fix on 

1 

the form u^u^^cj?. The vacuum oj mass is then m^^ = G^^cpo, where 0o denotes 
the vacuum glueball field (similarly ctq will be used for the sigma field). 
For the potential Vq we used [O 



Vg(0, a, tt) = In 41 - i - ^550^ In 



4 I 2 _2 
'0 / CTq 



(3) 



where ( = ^ and B and 6 are parameters. Here the logarithmic terms 
contribute to the trace anomaly: in addition to the standard contribution 
from the glueball field [|13|, |T^ there is also a contribution from the a field. 
Specifically the trace of the "improved" energy-momentum tensor is 

.. = 4M*.)-i:*.|| = 4.,„gy, (4) 

where $j runs over the scalar fields {</>, a, tt} and the vacuum energy, evac = 
— |i?0Q(l — 6). Guidance on the value of 6 is provided by the QCD trace 
anomaly which is proportional to the beta function. At the one loop level, 
with Nc colors and Uf fiavors, this is given by 

where the first number in parentheses arises from the (antiscreening) self- 
interaction of the gluons and the second, proportional tonj, is the (screening) 
contribution of quark pairs. This suggests a value of 6 = 4/33 for the present 
case with nj = 2 and Nc = 3 and we use this here, as in I, although nuclear 
properties are not sensitive to the precise value. The third term in the 
expression for Vq in Eq. (^) is needed to ensure that in the vacuum (p = (pQ, 
o" = (To and TT = 0, i.e. it provides for spontaneously broken chiral symmetry. 

2.1 Explicit Chiral Symmetry Breaking 

In order to generate a finite pion mass it is necessary to include explicit chiral 
symmetry breaking terms in the Lagrangian. To restrict the possibilities we 
recall SU{2) symmetry breaking at the quark level which is generated by 
VsB = rh{uu + dd), where m is an average mass. This gives a contribution 
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O^sb) to the trace of the energy- momentum tensor, since the operator is of 
dimension 3. We require our hadronic potential to have similar dimension. 
The double commutator of Vsb with the axial charge [Q^'', Vs^]] = 
'^SB^ab- We require an analogous relation to be satisfied in operator form for 
our hadronic potential. 

From Eq. the axial current and charge are 



It is straightforward to verify that there are three forms for the symmetry 
breaking which satisfy our criteria and we write the general potential 



In order to ensure that in the vacuum (j) = (j)Q, a = ao and tt = we need to 
add additional non-symmetry breaking terms to Vsb] we also subtract the 
vacuum value of these pieces. It is convenient to set e[ = eiaocpl, e'2 = €2(1^ 
and work with the ratio of the fields to their vacuum values, x = 0/0o and 
u = (j/aQ. Then 



Notice that V^b gives a (relatively small) contribution to the trace of the 
energy-momentum tensor. Specifically the contribution to the vacuum energy 
is 





- e2cr(o-2 + TT^) + e3NN . 



(7) 




e 



vac 




(9) 
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where for the latter equahty, involving the pion decay constant and mass, 
we have used the relations given in Subsec. 2.4 below. This is the standard 
result for the vacuum energy arising from the symmetry breaking term. 

Using the Lagrangian Ci — V^^, it is appropriate to give Lagrange's equa- 
tions for nuclei or nuclear matter at this point, since tt = in the mean field 



approximation. For the nucleon, the Dirac equation is of standard form [|15 
but now the effective mass is M* = ga + es = Mu + £3(1 — z/). The equations 
for the fields x{f^) = 0(^)/0O! = o"(^)/o"o and uJo{r) can be written 

(PlVx - 2[fio(2 -6) + e[ - Se'Jx + 2[Bo6 - 3e',]iy 

= 4Bo[xHlnx-Slnu) -X + H+ Bo6[x{i^' - x') + 2(x - u)] 
-mlulx - t\{2xv -X"- Xy'' + 2x) - Se'^fx' - X^'' - 2^ + 2v\ , 
alVv - {2BoS + e[ - 3e^)z/ + 2{BoS - Sea)^ 



gaops - Bq5 



A 



^-X^v + 2{v-x) 



x_ 

V 

-3e'2(z/2-xV + 2x-z/) , 
Pcjo - m>o = -gu^pB + ml{x^ - l)cJo + 4(^4)^^^^ , (10) 

where the densities ps = (NN) and Pb = (N'^N) can be expressed in terms 
of the components of the nucleon Dirac spinors in the usual way |T^. In 



eq. (|iy) we have made the definitions T) = + ^ J: and Bq = Bcf)^. The 
terms linear in the fields, i.e., the kinetic energy and mass terms, have been 
separated out on the left of these equations. In the case of nuclear matter 
the fields are constant and so the derivatives are zero. We note that the 
mass matrix is not diagonal and in order to solve the equations for nuclei 
it is necessary to go to a representation which is diagonal in the limit that 
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r oo, i.e. 1- This is discussed in the Appendix A. 

The energy-momentum tensor can be used to obtain the total energy of 
the system in the standard way Subtracting constants so that the energy 



is measured relative to the vacuum, we obtain 

Q^*^ POO f 

E = y] ea{2ja + 1) - 27r / dr r'^lgaoups + q^jUJoPb 

Jo I 

+2Bo[x\lnx-Slnu + i) - i] + lBo6[x\2u' - 3x') + 1] 
-mlx'ujl - 2{G,ru;', + le[[x'{x' - 2u + 2u') - 1] 

+ ie^(6z.V-3x'-2z.3-l)j. (11) 



In the first term on the right the ea are the Dirac single particle energies 
and ja is the total angular momentum of the single particle state. In nuclear 
matter this term becomes 4 (s'l^^o + \/k'^ + M*^^ and the integral in Eq. 
(p!Tf) is trivial since the fields are constant. 



Finally, for nuclei the photon field is included in the standard way [15 
and, since the contribution from the p field is small, it is handled in the 
simplest manner as in I. These fields are suppressed in the above equations. 

2.2 Axial Coupling Constant 

Eqs. (H) and (|) yield, in the vacuum, ga^ = M — 63 as the form of the 



Goldberger-Treiman relation |17, 21] which corresponds to a nucleon axial 
coupling constant qa = ^- It has been pointed out by Lee that this 
defect can be corrected by introducing an additional chiral invariant term 



in the Lagrangian. This has also been discussed by Akhmedov |T8[ in the 



context of the standard chiral model, which, as we have remarked, is unable 
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to reproduce the observed properties of finite nuclei. 

It is straightforward to motivate the additional term. With the definitions 

T, = a + iTT ■ T ; U = a + iir ■ T75 , (12) 

one observes that 

iV7^nt [d^u) N = iVK7^st [d^E) Nr + Nli.t. (a^st) Nl , (13) 

where Nl,r are the nucleon solutions of left- and right-handed chirality. Un- 
der global SUl{2) X SUr{2) transformations S ^ LER^ , Nl LNl and 
Nn RNr, so the expression (|13D is clearly invariant. 

Simply Hermitizing the expression (|13D yields N^f.d^' (u^u) N and, since 
the nucleon vector current is conserved, this can be written as a total diver- 
gence which gives no contribution to the action. A more interesting possibil- 
ity i is 



N 



= ^N-f^^T ■ [tt X a^TT + 75 (ad^Tz - Tzd^'a)] N , (14) 

where we have inserted a constant D and a factor of to ensure scale 
invariance. Since in a mean field calculation of nuclei or nuclear matter tt = 0, 
C2 will not modify the results obtained in I. It will, however, contribute to 
the axial vector coupling constant, qa, and the vrA^A^ coupling, gnNN, so that 
[y, |TB[ the Goldberger-Treiman relation q-^nnc^q = QaM — €3 is obtained. 



Before giving further details, the renormalization of the pion field needs to 
be addressed. 
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2.3 IT — ai Mixing 

The inclusion of the isotriplet vector and axial vector mesons (p and ai), along 
with the scalar and pseudoscalar mesons, in a linear chiral model results in 



significant complication. The details are discussed by Ko and Rudaz |]19 
(see also Gasiorowicz and Geffen [|10|) and their vector meson Lagrangian is 
given in Appendix B. When the vector mesons are included, the principle of 
vector meson dominance in the currents is implemented by requiring that, 
apart from explicit symmetry breaking terms, the Lagrangian be invariant 
under local SUl{2) x SUji{2) save for the p and Oi mass terms . This requires 
the replacement of the ordinary derivatives of the a and tt fields in Eq. (|^) 
by covariant derivatives: 

d^a A^a = d^a + /a^ ■ tt , 

d^TT A^TT = d^TZ + fp^XTZ- faRf, , (15) 

where / is the gauge coupling constant for the vector fields. This induces 
mixing between the vr and the Oi which results in a renormalization of the 
pion field. Considering just the kinetic energy terms for the tt and ai, the 
mass term for the ai and the mixing, we have 

Id^TT ■ d^Tz - ia^, • a'^'' + imf ■ a^ - fad^n ■ s!^ 

\z:d,7r ■ d'^Tz - la;, ■ a"^^ + >f a; ■ a"^ , (16) 

where a^j, = d^a^, — d^a^, the transformed field is 
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and 

Z: = l- ifa/mlf . (18) 

Here, and in the following, an asterisk indicates that the values of the a and 
fields in nuclear matter are used; thus m* denotes the mass of the ai in 
matter. Guided by the observation in I that nuclei prefer the p mass to be 
generated by coupling to the glueball field rather than the a field, we take 



the constants b and c of Ref. |]T9| to be equal so that the masses in matter 
are 

m:' = mf + {2c+l){far ; mf = mll^^J , (19) 

where the absence of an asterisk on the mass implies the vacuum value. Here 
we have implicitly assumed that the fields are replaced by their average values 
in nuclear matter, i.e. used the mean field approximation, in which case Z* 
is a renormalization and the physical pion field is defined by tt^ = y^Z*^7r. 

The redefinition (|l^ of the ai field in the Lagrangian (^31) of Appendix 
B and in the aiN coupling, —^fN'j^'j^TN ■ a^, leads to additional terms 
involving the a and tt fields. We can neglect those which are 0{tt'^) or 
0{dadTTY since they will not contribute here. The additional terms we need 
to consider are 

£3 = \{Zl-l) [d^TT ■ d^Tz + a-^N^^^T ■ Nd^iz] +L^d>'a7T-d^7r , (20) 



where we have included the renormalization of Eq. (|T^). Again, these terms 
will not affect nuclear matter at the mean field level. 
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We must also replace the derivatives in Eq. (JT^) by covariant derivatives 
yielding the modified form 



D 



(1 - Zl)a-^Tz{TZ ■ d^Tz] 



N . 



Thus our total Lagrangian is 



total 



(21) 



(22) 



2.4 Effective Couplings and the Pion Mass 

The axial current, A^, arises from the breaking of local SUl{2) x SUji{2) by 
the p and Oi mass terms and can be written [ff^] 



-cf ha-rz X 



a,. 



2(a^ ■ 77)77 - 77 



(23) 



Using this, we can determine the axial vector coupling constant in medium, 
g\, by a tree-level calculation of {N{k — q)\A^\N{k)) and identification of 
the coefficient of —N'~^^'~^^^tN in the limit g 0. This gives 



9a 



(24) 



Since in the vacuum < 1, the presence of the new Lagrangian £2 is 
essential if one is to obtain gA> ^■ 

Turning to the effective tiNN coupling —iglj^^N-^^T ■ tt^'N, where 77'' 
denotes the renormalized (physical) pion field, we find 



9nNN 



Ox Z* 



(25) 
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Here we have used the relation 



for on mass-shell nucleons. 

Using Eqs. (0) and (p3D, we have 



(26) 



(27) 



so that the effective pion decay constant in medium is /* = cr^Jz^ (we use the 
standard nomenclature here while noting that there would be additional con- 
tributions to the calculation of the decay of an actual pion in medium) . Then 
Eqs. 



and 



yield the Goldberger-Treiman [17| relation in medium 



(28) 



which is a consequence of the chiral symmetry employed. 

Finally we discuss the pion mass in symmetric (isospin zero) nuclear mat- 
ter. In addition to the contributions that are generated by straightforward 
application of Lagrange's equations to £totai of Eq. (|5), one must include 
[ pO|] the fluctuations in the axial density of nucleons that are induced by 
the presence of pions. This is represented by the propagator contribution in 
Fig. 1. We choose to evaluate the contribution to the pion mass at four- 
momentum g = (rather than identifying the pole of the propagator at 
q = 0) which yields g'^ps/{Z*M*). In this approximation the effective pion 
mass is given by 



*2 



(29) 
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In the vacuum this evidently becomes 

e'l + 6^ 



(30) 



We can now evaluate the pion-nucleon sigma term, which with our symmetry 
breaking is simply {N{p')\Vsb\N{p)), evaluated at t = — p'Y = 0. In the 
tree approximation, following the approach of Campbell , we obtain 



S(0) 



2gcroZ^ 



3m^7 + 



2e; 



using the definitions of Appendix A. 



(/^CTo - 70o) 



+ £3 



(31) 



3 Results 
3.1 Nuclei 

We shall examine the results obtained with and without the {uj^uj^Y term in 
the Lagrangian. For the latter case, the binding energy/particle (16 MeV) 
and density (p^* = 0.148 fm~^) of equilibrium nuclear matter are sufficient 
to fix the parameters and Bq for various assumed values of the symmetry 
breaking constants e'^, e'2 and £3. For the former, an additional parameter 
enters and we take G4,/ = 0.19. This is similar to the value used in I (where 
the glueball field was frozen, x = 1) and is judged to give the best results 
for the nuclear properties; smaller values of produce little effect, whereas 
significantly larger values preclude a fit to equilibrium nuclear matter. As in 
I, two solutions are found and we pick that which yields the most reasonable 
effective mass, M*^^, and compression modulus, K. The parameters are listed 
in Table 1 for the cases where 0^ = and 6*4 > 0. These parameters actually 
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correspond to our favored values for the symmetry breaking, with e[ chosen 
to give the pion mass (Eq. (pOD) for e'a = 0, and = —15 MeV. However 
the quantities hsted are insensitive to the values of these parameters since 
the scale of the symmetry breaking is small in comparison to the other scales 
involved. 

The first few quantities in Table 1 are similar to those discussed in I. 
Thus the vacuum energy is in rather good agreement with QCD sum rule 
estimates fH of (240 MeV)^, as in I, while the u coupling is comparable 



to that needed for the nucleon-nucleon interaction [Q. Since nuclei are not 
sensitive to the ratio ( = (po/cxQ we fix the value such that the larger scalar 



mass m> = 1.5 GeV in view of QCD sum rule estimates pj] for a dominantly 
glueball state. The lower mass m< ^ 0.5 GeV. The saturation effective mass 
ratio agrees with the recent QCD sum rule estimate [2^ of 0.73 ± 0.09. As 
discussed in I, the lower effective mass for G4 > improves the spin-orbit 
splittings in nuclei which are about 80% of the observed values, versus 65% 



when G4 = 0. Reference also determines the vector self-energy of a 
nucleon in nuclear matter, g^^ujo, to be 0.31 ± 0.06 GeV. We obtain 0.22 GeV 
with = and 0.28 GeV for 6*4 > which suggests that the latter yields 
the better phenomenology. As regards the compression modulus, K, current 



estimates are in the range 200-300 MeV ||2^ which again favors the case with 
6*4 > 0. Defining the third derivative of the binding energy per particle as 
^ ~ (^) (evaluated at equilibrium), we give in Table 1 the ratio S/K 

since the data indicate a linear relation between this quantity and K. 
For 6*4 > our point lies within the allowed error band, while for G4 = it 
is slightly off the band. 
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For nuclei we need to fix ctq and the values which give the best fit are 
listed in Table 1; a reduction in these values leads to poor agreement, but a 
modest increase does not substantially change the results presented in this 
section. The vacuum pion renormalization, Z^, and the gauge coupling, / 
can then be determined. Since the latter is the pNN coupling constant, it 



can be compared with the value of 5.3±0.3 given by Dumbrajs et al. 

our results are of reasonable magnitude. The remaining quantities in Table 

1 will be referred to below. 

Since the single particle levels and charge density distributions are quite 
similar to those obtained in I, we shall content ourselves with discussing the 
bulk properties of oxygen, calcium and lead. These are given in Table 2; 



here the appropriate corrections for cm. and finite size effects [|I5[ have 
been made. We also list the value of the sigma term from Eq. ([311) . With 
the choice (y'^^^/(47r) = 14.3 from Hohler the vacuum axial coupling 



contant, qa-, is determined and the values in Table 2 are seen to be quite 
reasonable. The parameter sets designated with a W in Table 2 correspond 
to the case with > 0. The cases labelled A correspond to no explicit 
symmetry breaking so, while qa is reasonable, the sigma term is necessarily 
zero; the remaining results are similar to those given in I. The other cases 
show the effect of various symmetry breaking parameters; here e'l and eg are 
constrained by Eq. (|30|) to give the vacuum pion mass. 

The symmetry breaking tends to decrease the binding energy and increase 
the radius or vice versa, but, as expected, the effects are fairly small. Simply 
choosing e[ ^ (case WB) improves the radii at the expense of some loss 
of binding. Also the sigma term becomes non-zero and is reasonable in 
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comparison to the estimate of Gasser et al. of 45 MeV with ~15% error. 
If we include e'2 (case WC), the radii are reduced and the sigma term becomes 
rather large, although the latter could be corrected by an appropriate choice 
of €3. We have also examined the choice = — ^e'l, but this significantly 
increase K and correspondingly substantially reduces the binding energies. 
We also note that with 63 7^ 0, the term in Vsb would dominate over the 
potential Vg as cr ^ 00, which would seem to be unphysical. We therefore 
favor setting €3 = 0. With this choice, we show in Table 2 the results obtained 
with various values of 63. A small negative value of ~ — 15 MeV appears to 
give the best results for nuclei and also for the sigma term. Set WF gives 
rather a good account of the radii with binding energies that are somewhat 
low, while set F improves the binding energies at the expense of the radii. 



3.2 ttN Scattering 



The Lorentz-invariant scattering matrix is conventionally written in the form 
Nip2)TbaN{pi), where 



Tba 

Here gi, a {q2, b) specify the four momentum and isospin component of the 
incoming (outgoing) meson. It is straightforward, though tedious, to obtain 
the amplitudes in terms of the usual Mandelstam variables s = (pi + gi)^, 
t = {qi — ^2)^ and u = {pi — 52)^- The plus amplitudes receive contributions 
from nucleon exchange in the s and u channels and from the exchange of 
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scalar mesons in the t channel. With on-mass-shell pions, we find 



AM 



'29 



+ (t - 27 + 



\lZ^{t-2^){t + ml) 



rj2 



M 



2 

9-knn 



1 1 



(33) 



Here m> and m< denote the scalar meson masses which are given in the 
Appendix A, along with the definitions for a,P and 7. Notice that is 
of standard form |20[| . 

For the minus amplitudes, nucleon exchanges in the s and u channels 
again contribute and in the t channel we need to consider p exchange. For 
the latter, the pTcir coupling follows from the Lagrangian quoted in Appendix 

as 



B and is given by Ko and Rudaz |T9 
C 



■'pTVTT 



2Z^mt 



2ml 



(34) 



where p^^ = df,p^ - d^p^. 
The pA^A^ vertex is 



pNN 



yNYr ■ p,N 



(35) 



In addition we have a four point NNtttt interaction in Eq. 
be considered. The resulting amplitudes are 



T!) which must 



A(-) = , 



B(-) 



2 

~9ttNN 



1 



p 



+ 



S-M2 



+ 



2^. 



9l 



2m2 



+ 



2m2 



(36) 



In terms of these amphtudes, the s-wave scattering lengths and effective 
ranges, and rQ^\ and the p-wave scatterings lengths, a^^, are EO 



(±) 



4M2 



- (2M + m^)S^^ 



(±) 



1 



2Mm. 



(M + m^) 



1 - 



A, 



(±) 



M + 



5, 



(±) 



(37) 



2M/ " \ 2M^ 
where Co = [d{A + 171^,3) /dt]t=o, -Do = [d{A + m^i?)/(9s]t=o, with s and t 
taken to be the independant variables. Also r]~^ = 47r(M + m.^)/M. The 
subscript implies that the amplitudes are to be evaluated at threshold:- 
t = 0, s = (M + m^)2 and u = {M - m^f. 

In order to calculate the scattering amplitudes in the most general case, 
additional parameters are needed. The parameter D is determined from 
Eqs. ( plD and (^5]) in the vacuum and is listed in Table 1. The quantity c in 
Table 1 is deduced from the value of and Eqs. (0) and (p!9D in vacuum. 
The parameter Kg listed in Table 1 is obtained by fitting the decay width, 
r(p tttt) = 151.2 ± 1.2 MeV H], using Eq. (||) with the p on the mass 
shell. The effective pim coupling constant of 6.05 is of similar magnitude 
to the pNN coupling / of Table 1, which is in accord with the notion of 
p universality. Since significant cancellations may occur in calculating the 
scattering lengths and effective ranges, Eq. (|37D, we prefer to deduce Aq, -Bq, 
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Co and Do from the experimental data and these are hsted in the first row 
of Table 3. The A resonance contributions to the amplitudes, constructed 
such that no double counting occurs, have been given by Hohler p^. We list 



these in Table 3, along with the difference which we would hope to reproduce 
with our model. (For definiteness we note that the off-shell parameter Z 
appearing in the vrA^A interaction was chosen to be |.) We first give results 
obtained with our simplest Lagrangian, Ci. Since gnNN = in this model, 
the results are quite poor. If we also include £2, this allows us to fit gj^NN by 
adjusting D (here is still 1) and this gives much better agreement with the 
data. The remaining cases in Table 3 correspond to the full model and the 
notation follows that of Table 2. With no explicit symmetry breaking (cases 
A and WA), the sum of aI^^ and itlt^B^q^^ must cancel so that there is no 
contribution to Oq'^ in the chiral limit. In fact the sum should be negative, 
as it is for most of the cases displayed. The sensitivity to the symmetry 
breaking is not great, but there appears to be a weak preference for cases F 
and WF and we recall that these were also favored from the results of Table 
1. 

The question arises as to how well the scattering parameters test our 
model. In this connection it is useful to compare with the results of Matsui 



and Serot [|16], ^ who used both the standard chiral model and QHD-II. In 
all models Qt^nn is taken from experiment so that the predictions for B^^^ 
and D^^^ are the same. Also Aq^'' will only differ from —771^^3^^ because of 



symmetry breaking terms, which are necessarily small. In Ref. [16, Cg^^ 



is 2.0 (3.0) in the chiral (QHD-II) model where the scalar mass is 770 (520) 
MeV. Thus one distinct feature of the present model is that we get about the 
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right value for Cq with m< ~ 500 MeV. As regards the (— ) amphtudes, 
Dq~^ is fixed once gjrNN is given and to a large extent Co""* is also; thus 
all models give much the same result. However in QHD-II, where the p is 
included, 171.^3^ ^ = 3.1, whereas our value is smaller, due to the constant 
term and the more complicated form for the p contribution, which brings us 
closer to the desired value. 



3.3 Density Dependence of the Parameters 

In our Lagrangian there are two scales, x ^i^nd uj which give the ratio of 
the glueball and sigma fields to their vacuum values. Their dependence on 
density is shown in Fig. 2 for the parameters of Table 1. The glueball field 
changes very little which implies that at saturation the uj and p masses are 



reduced by only 3%. By contrast QCD sum rule estimates suggest a 
reduction of roughly 20 ± 10%. Brown and Rho have also argued for 
a similar scaling. An effect of this magnitude would suggest scaling these 
vector masses with z/ instead of x- Then, however, it is not possible |^ ^ || 
to produce acceptable fits to the properties of nuclei, for instance the charge 
density distributions start to develop oscillations in disagreement with the 
data and these probably signal ||3^ the onset of a situation where the ground 
state of nuclear matter is no longer uniform. 



The pion decay constant in medium is given by f*/ = lyy Z*/Zt,. The 
behavior of the pion renormalization, Z*/Z.,r, is shown in Fig. 3. The change 
in Z* in going from the vacuum to saturation is a 7-11% effect, depending 
on the parameters chosen, so that f*/ — M*/M as suggested by Brown 
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and Rho [Q. Specifically the pion decay constant ratio at saturation is 0.74 
(G'4 = 0) or 0.66 ((74 > 0) which compare quite well with the effective masses 
given in Table 1. A similar level of agreement is obtained if one also includes 



a factor of Qa- The ratio of the axial coupling constants, given in Fig. 
4, shows the expected reduction with density. The data indicate [^] that 
g\ in equilibrium nuclear matter should be close to 1, which is a little lower 
than our value of 1.1. 

The quantities discussed thus far are insensitive to the value of the explicit 
symmetry breaking parameters. This is not the case for the effective pion 
mass, however, as we show in Fig. 5. At low densities, where the linear 



approximation is sufficient, the behavior is given by the expression ||3^ , p8 
m*/m^ = 1 — ^Tia'o^ pBlivTLT^vTLji), where m/j is the reduced mass of the 
pion-nucleon system. (This can also be obtained in leading order from Eqs. 
(|29D, ([33| ) and (|37D .) The mass decreases for set WD since ag"*"^ is positive, 
whereas is should be negative leading to an increased pion mass in matter. 
The pionic atom data indicate a 10% increase at saturation and this 
suggests that €3 should be zero or small and negative in agreement with the 
previous indications. 



4 Conclusions 

The effective Lagrangian of I, which embodied broken scale and sponta- 
neously broken chiral symmetry as suggested by QCD, gave a good account 
of the properties of finite nuclei as well as nuclear matter. We have extended 
this Lagrangian by allowing for the explicit breaking of chiral symmetry and 
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by including an additional term which permits us to obtain an axial coupling 
constant, qa-, which is close to the physical value of 1.26 in vacuum. We have 
also adjoined this Lagrangian to a chiral Lagrangian for the p and ai vector 
mesons while maintaining reasonable properties for these particles. This 
required a renormalization of the pion field. The final Lagrangian yielded 
the Goldberger-Treiman relation at finite density, as well as in the vacuum. 

The effect of explicit chiral symmetry breaking on nuclear matter and 
nuclei was, for the most part, small so that their properties continued to 
be reasonably well described. As in I, the presence of a term {uo^uj^Y 
the Lagrangian improved the phenomenology. The low energy vrA^ scattering 
parameters were found to be in good agreement with the data once the effect 
of the A resonance was removed, since this is not included in our model. As 
regards the modifications of the various quantities in going from the vacuum 
to equilibrium nuclear matter density, the substantial reduction in the p and 
uj masses predicted by QCD sum rules remains puzzling since the reduction 
we obtain is very slight. Further it appears that a more substantial reduction 
yields poor results for nuclei P, ^, On the other hand, we find that the 
axial coupling is quenched to a value close to 1, as indicated experimentally, 
and the increase in the pion mass suggested by pionic atom data is obtained 
if the constant in the symmetry breaking potential e^NN is zero or very 
small and negative. The latter was favored by other data, particularly the 
pion-nucleon sigma term. Thus we conclude that our Lagrangian, in which 
chiral symmetry is realized in a linear fashion, embodies the essential physics 
necessary to describe quite a wide variety of low energy data. 

We acknowledge partial support from the Department of Energy under 
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Appendix A: Scalar Field Mass Matrix 

Setting (T = (To — o" and (p = (j)o — (f), we denote the 2x2 mass-squared" matrix 
for the fluctuating parts of the scalar fields, a and 0, by 

2a -2(3' 
-2(3 27 

The matrix elements can be read off from the field equations ([T0|) : 



(38) 



_ 2^0,5 + e[ - 3e^ _ _ Bpd - 3e^ _ _ Bo{2 - 6) + e[ - 3e'^ 



V^Ot^O V^O 



Defining r = ^{'-) — + 4/^2, the eigenvalues and eigenvectors are 
m> = a + 7 + r ; = aa — b(j) 

m\ = a + 7 — r ; = ba + acj) , (40) 

where a/6 = (a - 7 + r)/{2f3) = 2(3 /{-f - a + r) and a'^ + = 1. 

Appendix B: The Lagrangian for the Vector 
Fields 



Here we quote, for reference, the Lagrangian of Ko and Rudaz [|r^ for the 
vector and axial vector mesons, p and ai. They define left and right gauge 
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fields 

= K/^m + ^m) ■'^ ' '^/^ = •'^ ' (41) 

and field strength tensor 

If^v = d^l„ - dj^ - if[l^, ly] , (42) 

with an analogous definition for r^^. The Lagrangian is then 

Sector = -\ Tr[/2^ + r^J + \ml Tr[Zj + rj] 

+ \hp TrpSt] Tr[/J + rj] - cf Tr[/^Sr^St] 

-i^ Tr[/^,A^'SA''St + r^.A^STA'-S] , (43) 

using the definitions of Eqs. (|12D and (pISl). For simplicity we focus here 



on the minimal model of Ref. |19] which neglects a further term in the 
Lagrangian ((^e = 0). In the present work we have set 6 = c (see Table 1 
for the values) and generated the dimensionful constants in a scale invariant 
way, i.e. replaced ttiq by m^pi^cf) / (p^Y and i^Q/m?^ by (fi;6/m^)(0o/0)^. 

We should briefly comment on the vector meson phenomenology with 
our chosen parameters. The predictions for the widths and the pion radius 
discussed in Ref. [|19| are given in Table 4 for our parameters (note that the 



corrected formula [pQ] , ^ is used for the ai pn decay width). In general 
the results are reasonable, although the ai —>■ tt'-/ width is too large by a 
factor of two. This width is sensitive to the presence of a (q term in the 
Lagrangian which leads to a renormalization of the p and ai fields. However, 
inclusion of this term to improve the phenomenology is beyond the scope of 
the present work. 
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Table 1 

Values of the parameters and derived quantities in nuclear 



Quantity 


G4 = U 


(j4/fi'c^ = 0.19 


1, 11/4 /'l\/r„"^;"^ 


ZOD 


zzl 




10. 


lz.4 


C = 0O/C^O 


1.3 


1.4 




0.70 


O.o4 


K (MeV) 


388 


317 


S/K 


6.6 


5.4 


(To (MeV) 


110 


102 




0.71 


0.83 


f 


5.97 


4.95 


D 


1.36 


1.10 


c 


0.57 


1.30 


Ke 


-1.52 


-1.90 
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Table 2 

Bulk properties of nuclei, sigma term and vacuum axial vector coupling 
constant with various explicit symmetry breaking parameters. 







4 


£3 


9A 


E(0) 




BE 
A 


' en 


Ca 

BE 
A 


' en 


Pb 

BE 
A 


' on 




(MeV) 




(MeV) 




(MeV) 


(MeV) 


fm 


(MeV) 


fm 


(MeV) 


fm 




Expt. 


1.26 


7.98 


2.73 


8.55 


3.48 


7.86 


5.50 


A 











1.33 





7.35 


2.64 


7.96 


3.41 


7.30 


5.49 


F 


138 





-15 


1.31 


41 


7.40 


2.63 


7.99 


3.41 


7.34 


5.49 


WA 











1.33 





7.02 


2.68 


7.75 


3.45 


7.33 


5.51 


WB 


138 








1.33 


71 


6.58 


2.72 


7.43 


3.47 


7.18 


5.52 


wc 


138 







1.33 


101 


7.46 


2.65 


8.10 


3.42 


7.51 


5.51 


WD 


138 





50 


1.38 


120 


6.36 


2.73 


7.26 


3.48 


7.06 


5.53 


WE 


138 





-50 


1.27 


24 


6.43 


2.73 


7.31 


3.48 


7.13 


5.52 


WF 


138 





-15 


1.31 


56 


6.58 


2.72 


7.45 


3.47 


7.20 


5.52 
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Table 3 

Low energy pion scattering parameters. 





^0 


m b!\^^ 




^0 


^0 












V'^TT ) 


V"'Tr J 


( m^'^] 

V"'7r J 


\"''1T J 




\"''1T J 


( m~'^] 

V"'7r J 


F,ynt 




— S9 5+n 4 


9 88+n 02 


2 65+0 04 


— 1 7+0 4 


1 2 0+0 4 


—1 76+0 02 


— 1 1 3+0 04 


A 

i— i I'-'M 


7 n 


— 5 4 


n Q8 


65 


— 10 7 

X W . 1 


1 s 

X w . o 


—0 64 


—0 85 


Til TT PTPTl PP 




—27 1 


1 8Q 


2 00 


1 

W . X 


1 7 

X . ( 


X . X X 


—0 28 


''-1 


ID 7 




1 18 

X . X (J 


79 




8 


—0 39 


—0.12 


A + /:2 


26.4 


-26.4 


1.76 


1.97 




0.8 


-0.97 


-0.29 


Case A 


26.4 


-26.4 


1.64 


1.97 




1.1 


-0.97 


-0.29 


Case F 


25.9 


-26.5 


1.79 


1.97 




1.2 


-1.09 


-0.29 


Case WA 


26.4 


-26.4 


1.87 


1.97 




1.1 


-0.97 


-0.28 


Case WB 


25.8 


-26.5 


2.01 


1.97 




1.1 


-1.10 


-0.29 


Case WC 


26.4 


-26.5 


1.99 


1.97 




1.1 


-1.10 


-0.29 


Case WD 


26.6 


-26.5 


2.00 


1.97 




1.0 


-1.10 


-0.29 


Case WE 


24.9 


-26.5 


2.06 


1.97 




1.3 


-1.10 


-0.29 


Case WF 


25.6 


-26.5 


2.02 


1.97 




1.2 


-1.10 


-0.29 
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Table 4 

Predicitions for the vector mesons with our parameters. 



Quantity 


= ^4/^,^ = 0.19 


Experiment 


r(ai ^ ptt) (MeV) 
ai ^ pn D / S ratio 
r(ai ^ 7r7) (MeV) 
r(pO ^ e+e") (keV) 
{r')l (fm) 


306 601 
-0.25 -0.17 
1.10 1.19 
4.82 7.00 
0.63 0.69 


~ 400 
-0.14 ±0.03 41| 
0.64 ±0.25 
6.77 ±0.32 
0.66 ±0.01 



Figure Captions 

Figure 1. Nucleon loop contribution to the pion propagator. 

Figure 2. Ratio of glueball and sigma fields to their vacuum values, x and 

z/, as a function of density. The parameters F and WF are specified in the 

text. Equilibrium nuclear matter density is denoted by p^*. 

Figure 3. Ratio of the pion renormalization Z* in matter to the vacuum 

value Zt, clS db function of density. 

Figure 4. Ratio of the axial coupling constant g\ in matter to the vacuum 
value qa as a function of density. 

Figure 5. Ratio of the pion mass m* in matter to the vacuum value 
function of density. The parameter sets are specified in the text. 
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